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P U  RPOS E :  The  Coastal  and  Hydraulics  Engineering  T echnical  Note  (CHETN)  described  herein 
provides  a  new  formula  for  estimation  of  irregular  wave  runup  on  rough,  impermeable  slopes.  The 
runup  guidance  is  based  on  the  wave  momentum  flux  parameter  described  in  CHETN-III-67 
(Hughes  2003a).  Sample  calculations  illustrate  application  of  the  formula. 

BAC  KGROU  N  D :  Estimates  of  maximum  wave  runup  on  rough,  impermeable  sloping  structures 
(riprap  revetments)  are  necessary  to  determine  whether  overtopping  will  occur  for  a  specified  wave 
condition  and  water  level.  Design  fonnulas  were  originally  developed  based  on  theory  and  small- 
scale  laboratory  experiments  using  regular  waves.  As  laboratories  acquired  the  capability  to  generate 
more  realistic  irregular  waves,  improved  wave  runup  formulas  were  proposed  based  on  wave 
parameters  representative  of  the  irregular  wave  train.  However,  unlike  regular  waves  that  result  in  a 
single  value  of  maximum  wave  runup,  irregular  waves  produce  a  runup  distribution.  Thus,  it  was 
necessary  for  the  runup  formulas  to  detennine  a  representative  parameter  of  the  wave  runup 
distribution.  Presently,  the  most  common  irregular  wave  runup  parameter  is  Ru 2%.  This  parameter  is 
defined  as  the  vertical  distance  between  the  still-water  level  (swl)  and  the  elevation  exceeded  by 
2  percent  of  the  runup  values  in  the  distribution.  In  other  words,  for  every  100  waves  running  up  a 
slope,  two  waves  would  have  a  runup  elevation  exceeding  the  level  estimated  by  Ru 2%. 

Irregular  wave  runup  design  guidance  for  rough,  impermeable  slopes  given  in  the  Coastal 
Engineering  Manual  is  based  on  runup  measurements  acquired  during  irregular  wave  rock  armor 
stability  experiments  conducted  at  Delft  Hydraulics  and  reported  by  van  der  Meer  and  Stam  ( 1 992). 
The  guidance  is  presented  as  two  empirical  fonnulas  with  each  formula  valid  over  a  specific  range  of 
Iribarren  numbers.  The  formulas  in  the  Coastal  Engineering  Manual  for  the  2-percent  runup 
elevation  are  as  follows: 


for  1.0  <  £,om  <1.5 
for  1.5  <  E,om 


(1) 


with  the  Iribarren  number  defined  as: 


tan  a 


(2) 


'om 
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and 


R-u2%  =  vertical  runup  distance  exceeded  by  2  percent  of  wave  runups 
Hmo  =  zeroth-moment  energy-based  significant  wave  height 
corn  =  deepwater  Iribarren  number  based  on  mean  period  Tm 
Lom  =  deepwater  wavelength  [=  (gfln)  Tm2] 
g  =  gravitational  acceleration 

Tm  =  mean  wave  period  associated  with  wave  spectrum 
tan  a  =  structure  slope 


The  mean  wave  period  is  used  instead  of  the  peak  spectral  wave  period,  Tp,  in  the  runup  formulas  to 
accommodate  different  widths  of  the  wave  spectrum.  However,  in  some  cases,  design  wave 
conditions  are  specified  in  tenns  of  Tp,  so  it  is  necessary  to  give  an  estimate  of  Tm.  The  ratio  of  TpITm 
varies  for  different  spectrum  types  as  follows: 

Tp  _  f 1  ’ 1 5  - 1 21  for  J0NSWAP  spectra 

Tm  [1 .22  - 1 .4 1  for  Pierson  -  Moskowitz  spectra 

van  der  Meer  and  Stam  (1992)  gave  the  range  Tp  lTm  =  cop  / com  =  1.1  -  1.2  when  they  converted 
Hunt’s  runup  equation  to  use  an  Iribarren  number  based  on  peak  spectral  period,  cop,  instead  of  com- 
Using  the  median  value  in  this  range  of  cop  / com  =  1.15,  Equation  2  can  be  converted  to  a  form  that 
uses  Tp  rather  than  Tm,  i.e., 

R  f 0.835  £ 

K„2%  =  I  h°P 

Hmo  1 1.10  (£,op)0'46 


for  1.15  <£*<1.72 

for  1-72  <£op 


where 


tan  a 


(5) 


and  Lop  is  the  deepwater  wavelength  based  on  peak  period,  Tp. 

Figure  1  presents  the  original  data  of  van  der  Meer  and  Stam  (1992)  along  with  additional  laboratory 
observations  reported  by  Ahrens  and  Heimbaugh  (1988).  The  observed  values  of  dimensionless 
runup  ( RuiWHmo )  were  plotted  versus  Iribarren  number  based  on  peak  period  ( qop)  for  structure 
slopes  of  1:2,  1:3,  and  1:4.  The  experiments  of  van  der  Meer  and  Stam  were  for  structures  with 
relatively  deep  water  at  the  toe,  and  the  maximum  value  of  relative  wave  height  was  about 
Hmo/h  =  0.25.  Ahrens  and  Heimbaugh’ s  experiments  used  shallower  water  depths  with  the  maximum 
value  of  relative  wave  height  of  about  Hmo/h  =  0.64.  The  solid  lines  in  Figure  1  were  plotted  using 
Equation  4.  Generally,  the  empirical  runup  equations  represent  the  data  well. 
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Van  der  Meer  and  Stam  (1992)  2%-Runup  Equations 


Iribarren  Number  - £ 

^op 


Figure  1 .  Irregular  wave  runup  guidance  given  in  the  Coastal  Engineering  Manual  for  rough, 
impermeable  slopes 


NewWave  Runup  Equation:  Hughes  (2003a,  2004a)  presented  anew  nondimensional  param¬ 
eter  representing  the  maximum  depth-integrated  wave  momentum  flux  that  occurs  in  progressive 
water  waves.  The  parameter,  referred  to  as  the  wave  momentum  flux  parameter  ( Pmf ),  was  defined 
as 


MF 


r  MF  a 

pgh2 


max 


where 

Mf  =  depth-integrated  wave  momentum  flux 
p  =  fluid  mass  density 
g  =  gravitational  acceleration 
h  =  water  depth 

Because  (MF )nYdX  has  units  of  force  per  unit  wave  crest  length,  it  was  argued  that  maximum  depth- 
integrated  wave  momentum  flux  would  provide  a  good  characterization  of  wave  processes  at  coastal 
structures. 
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Hughes  (2003a,  2004a)  established  an  empirical  equation  for  estimating  the  wave  momentum  flux 
parameter  for  finite  amplitude,  nonlinear  waves  based  on  a  numerical  solution  technique  (Fourier 
approximation).  The  resulting,  purely  empirical  equation,  was  given  as 


(  Mf  ^ 

v  P 'gti  J 


Ar 


(  h  XA' 


\§r  J 


where 


A0  =  0.639 


A1  =0.180 


H 

\h  j 

H 


V  h  j 


(6) 


(V) 

(8) 


and  H  and  T  are  the  regular  wave  height  and  period,  respectively.  More  infonnation  and  a  sample 
calculation  of  wave  momentum  flux  are  given  in  CHETN-III-67  (Hughes  2003a). 

Hughes  (2003b,  2004b)  derived  a  simple  generic  expression  for  wave  runup  in  tenns  of  wave 
momentum  flux.  He  argued  that  the  weight  of  fluid  within  the  hatched  wedge  area  (ABC)  of 
Figure  2  is  proportional  to  the  maximum  depth-integrated  wave  momentum  flux  of  the  wave  just 
before  it  reached  the  toe  of  the  structure  slope.  The  resulting  equation  was  in  the  form 


R 

h 


CF{  a) 


Mf 

pgh1 


v  1/2 


(9) 


where  R  is  the  maximum  vertical  runup  from  swl,  and  C  -F(a)  is  an  unknown  function  of  structure 
slope  angle.  For  convenience  the  subscript  “max”  has  been  dropped  from  the  wave  momentum  flux 
parameter. 

In  the  Equation  9  runup  equation,  relative  runup  (R/h)  is  directly  proportional  to  the  square  root  of 
the  wave  momentum  flux  parameter.  (Note  that  representing  the  runup  sea  surface  slope  as  a  straight 
line  is  an  approximation  and  may  not  be  fully  appropriate  for  nonbreaking  waves  on  steep  slopes 
where  the  water  surface  has  pronounced  curvature.) 

Applying  Equation  9  to  irregular  wave  runup  requires  that  regular  wave  height  and  period  (H  and  T) 
used  to  estimate  the  wave  momentum  flux  parameter  using  Equations  6,  7,  and  8  be  replaced  with 
representative  irregular  wave  parameters  ( Hmo  and  Tp).  This  substitution  does  not  imply  that  an 
equivalence  exists  between  values  of  wave  momentum  flux  parameter  calculated  for  regular  and 
irregular  waves;  it  only  provides  a  convenient  standard  for  application  with  irregular  waves  when 
establishing  empirical  relationships.  For  irregular  waves,  always  check  that  the  depth-limited  wave 
height  does  not  exceed  the  rule-of-thumb  Hmo  <  0.6  h. 
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Summary  of  Irregular  Wave  Runup  Prediction  for  Smooth  Slopes:  Hughes  (2003b) 
recommended  the  following  equations  for  prediction  of  irregular  wave  runup  on  smooth, 
impermeable  slopes  for  nonbreaking  (surging/collapsing)  waves  and  breaking  (plunging/spilling) 
waves  as  delineated  by  values  of  Hmo  ILP. 

Nonbreaking  (surging/collapsing)  waves  ( HmoILp  <  0.0225): 


^»2% 

h 


1.75 


l-gd1-3™101) 


Mp 

?gh  J 


\  1/2 


for  1.0  <cot  a  <4.0 


(10) 


Breaking  (plunging/spilling)  waves  (Hmo/Lp  >  0.0225) 

R 


i/2% 

h 


=  1.75 


l  +  e-(°-47cota) 

f  mf  \ 

L  J 

Ipg^2  J 

1/2 


for  1.5  <  cot  a  < 4.0 


(11) 


The  formulas  are  for  structure  slopes  steeper  than  1 :4  (cot  a  =  4). 

More  recently  Hughes  (2004b)  modified  Equation  1 1  for  breaking  (plunging/spilling)  waves  to 
include  irregular  runup  on  smooth,  impermeable  slopes  as  mild  as  1:30.  The  newer  equation  was 
given  as: 
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Newer  equation  for  breaking  (plunging/spilling)  waves  (H„JLP  >  0.0225): 


^»2% 

h 


4.4  (tana)07 


s  \  1/2 

A  mf  ^ 


pgh2 


for  1.5  <  cot  a  <  30 


(12) 


Equation  12  covers  a  much  broader  range  of  slopes,  but  it  is  slightly  less  accurate  than  Equation  1 1 
when  applied  to  structure  slopes  in  the  range  1.5  <  cot  a  <  4.0.  Nevertheless,  having  one  equation 
that  predicts  irregular  wave  runup  over  such  a  wide  range  of  slopes  is  convenient,  and  it  supports  the 
simple  concept  used  to  derive  the  runup  equation. 

Irregular  Wave  Runup  Prediction  for  Rough,  Impermeable  Slopes:  Slope  roughness  will 
reduce  the  2-percent  runup  level  predicted  using  the  equations  for  smooth,  impermeable  slopes 
(Equations  10  and  11).  One  engineering  approach  is  to  multiply  the  runup  estimates  for  smooth, 
impermeable  slopes  by  a  reduction  factor  to  account  for  various  types  of  slope  roughness.  The 
Coastal  Engineering  Manual  contains  reduction  factors  as  summarized  by  de  Waal  and 
vanderMeer  (1992)  based  on  Dutch  experience,  and  the  reduction  factor  for  rock  and  riprap 
structures  ( 1  and  2  layers)  varies  between  0.5  and  0.6.  Rock  and  riprap  structures  impede  runup  not 
only  by  slope  roughness,  but  also  by  penneability  of  the  riprap  and  any  underlayers  placed  over  the 
impermeable  slope.  Voids  are  constantly  filling  and  draining  with  the  wave  runup/rundown  cycle, 
and  the  effect  of  riprap  permeability  will  vary  with  wave  period.  Thus,  using  a  single  constant  to 
represent  runup  reduction  associated  with  rough,  impermeable  riprap  slopes  is  simplistic.  However, 
this  approach  is  justified  by  the  success  of  past  runup  estimation  fonnulas. 

The  original  riprap  slope  runup  data  of  van  der  Meer  and  Stam  (1992)  and  Ahrens  and  Heimbaugh 
(1988)  were  plotted  as  a  function  of  the  wave  momentum  flux  parameter  to  ascertain  any  difference 
between  breaking  (plunging  and  spilling)  waves  on  the  slope  and  nonbreaking  (surging  and 
collapsing)  waves.  No  strong  trend  was  evident  for  rough  slopes,  whereas  a  distinct  difference  was 
apparent  for  runup  on  smooth  slopes  (Hughes  2003b,  2004b).  Because  the  majority  of  the  laboratory 
data  represented  waves  that  break  on  the  slopes,  the  newer  runup  equation  (Equation  12)  for 
plunging/spilling  waves  was  used  to  estimate  an  appropriate  reduction  factor  for  runup  on  rough, 
impermeable  slopes.  This  method  is  appealing  because  it  retains  the  conceptual  runup  model  shown 
in  Figure  2.  The  best-fit  reduction  factor  representing  the  mean  of  the  data  was  0.505,  which  is  at  the 
lower  end  of  the  reduction  factor  range  for  rock-annored  slopes  given  by  de  Waal  and  van  der  Meer 
(1992). 

Applying  the  reduction  factor  to  Equation  12  results  in  the  following  new  equation  for  estimating 
irregular  wave  runup  on  rough  (riprap)  impermeable  slopes  for  both  breaking  and  nonbreaking 
incident  wave  conditions. 


New  equation  for  irregular  wave  runup  on  rough,  impermeable  slopes: 


^»2% 

h 


4.4  (tana)07 


Mf 

pgh2 


,  1/2 


(0.505) 


for  2.0  <  cot  a  <  4.0 


(13) 


Figure  3  compares  estimated  values  of  relative  runup  ( RuivJh )  to  measured  values.  Reasonable 
correspondence  was  expected  because  these  were  the  same  data  used  to  determine  the  appropriate 
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reduction  factor.  The  data  of  van  der  Meer  and  Stam  (1992)  are  clustered  at  low  values  of  Ruia/Jh 
because  the  experiments  were  conducted  for  structures  in  relatively  deep  water.  In  deeper  water, 
incident  wave  nonlinearity  is  less,  and  the  waves  bear  more  resemblance  to  sinusoids. 

A  direct  comparison  between  the  new  runup  formula  (Equation  13)  and  formulas  given  in  the 
Coastal  Engineering  Manual  (Equation  4)  for  the  two  laboratory  data  sets  is  made  in  Figures  4  and 
5.  The  new  runup  equation  based  on  wave  momentum  flux  (Figure  4)  gave  a  reasonable  estimate  of 
the  measured  2-percent  runup,  but  the  van  der  Meer  and  Stam  (1992)  formulas,  converted  to  use 
peak  period,  Tp,  gave  slightly  better  results  (Figure  5).  Most  of  the  difference  between  the  two 
methods  stems  from  data  obtained  for  the  1 :2  rough  slope,  particularly  the  data  of  van  der  Meer  and 
Stam  (1992)  that  appear  to  be  entirely  nonbreaking  waves  with  all  values  of  %op  >  4.4. 


Predicted  vs.  Measured  R  _/  h  -  Momentum  Flux 

112% 


Estimated  R  /  h 

112% 


Figure  3.  Predicted  versus  measured  Ru2%/b  for  rough,  impermeable  plane 
slopes 
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Estimated  Runup  -  Momentum  Flux 


Estimated  Ru2%  -  (m) 


Figure  4.  Runup  estimation  using  wave 

formulas  momentum  flux  formula  (Eq.  13) 


Estimated  Runup  -  Van  der  Meer  &  Slam  1992 


Estimated  Ru2%-(m) 


Figure  5.  Runup  estimation  using  Coastal 

Engineering  Manual  converted  to  use 
peak  period,  Tp  (Eq.  4) 


van  der  Meer  and  Stam  ( 1 992)  determined  their  empirical  runup  equation  by  a  direct  curve  fit  to  the 
data.  This  procedure  was  able  to  represent  the  nonbreaking  wave  data  better  than  the  method  of 
applying  a  uniform  reduction  factor  to  the  estimates  for  the  corresponding  runup  on  smooth  slopes 
using  the  plunging/ spilling  runup  equation  based  on  wave  momentum  flux.  The  nonbreaking  wave 
runup  data  were  isolated  and  fitted  with  a  separate  curve,  and  this  improved  the  wave  momentum 
flux  estimates  to  equal  those  of  van  der  Meer  and  Stam  (1992).  However,  the  improvement  was  not 
significant  enough  to  justify  having  two  formulas  when  a  single  formula  gives  estimates  comparable 
to  those  of  van  der  Meer  and  Stam. 


EXAMPLE:  IRREGULAR  WAVE  RUNUP  ON  ROUGH,  IMPERMEABLE  SLOPES: 


Find:  The  vertical  runup  distance  from  the  swl  which  is  exceeded  by  only  2  percent  of  the  waves 
(i.e.,  Ru 2%)  for  structure  slopes  of  1:2  and  1:4  (tan  a  =  0.5  and  0.25). 

Given: 


h  =  20  ft1 
Tp  =  9.0  sec 
Hmo  =  8  ft 
g  =  32.2  ft/sec2  - 
tan  a  =  0.5,  0.25  - 


water  depth  at  toe  of  slope 

wave  period  associated  with  spectral  peak 

zeroth-moment  significant  wave  height 

gravitational  acceleration 

structure  slope 


1  To  convert  feet  to  meters,  multiply  number  of  feet  by  0.3048. 
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Calculate  wave  momentum  flux  parameter:  First  calculate  values  of  relative  wave  height  and 
relative  depth  as 


H  _  8  ft 
T  ~  20  ft 


gT2  (32.2  ft/sec2)(9  sec)2 


=  0.0077 


Next,  find  the  values  of  the  coefficient  A0  and  A  \  from  Equations  7  and  8,  respectively,  i.e., 


A0  =  0.639  —  =  0.639  (0.4)^  =  0.0998 

\  h  J 


4  =0.180  —  =  0.180  (0.4)  =  0.2576 


Finally,  the  nondimensional  wave  momentum  flux  parameter  is  calculated  from  Equation  6  as 


J—  =  a 

,  2  ^0  rji2 


=  0.0998  (0.0077)  °-2576  =  0.35 


Calculate  runup  for  1:2  slope: 

For  rubble-armored  impermeable  slopes,  the  appropriate  runup  equation  for  both  surging/collapsing 
and  plunging/spilling  breaker  types  is  given  by  Equation  13  as 


^^  =  4.4  (tanaf7  I  ]  (0.505) 

h  K  >  [p  gh2 


First  check  that  the  structure  slope  falls  within  the  range  of  applicability  for  Equation  13,  i.e., 

cot  a  =  — - —  =  — —  =  2.0 
tan  a  0.5 

which  is  within  the  range  of  2.0  <  cot  a  <  4.0. 

The  nondimensional  relative  2-percent  runup  is  found  as: 


^^  =  4.4  (tana)07  [  ]  (0.505)  =  4.4  (0.5)°'7(0.35)1/2(0.505)  =  0.81 

h  Ip  gh2 


and  the  dimensional  2-percent  runup  is 

Ru2%  =0.81/7  =  0.81  (20  ft)  =  16.2  ft 
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For  comparison,  the  present  Coastal  Engineering  Manual  method  given  by  Equation  4  estimates 
runup  to  be  Rui%  =  15.9  ft. 


Calculate  runup  for  1:4  slope: 

After  checking  that  the  1 :4  structure  slope  falls  within  the  range  of  applicability  for  the  plunging/ 
spilling  wave  runup  equation,  the  nondimensional  relative  2-percent  runup  is  found  from 
Equation  13  as 


^»2% 

h 


4.4  (tana)07 


Mf 

pgh2 


,1/2 


(0.505)  =  4.4  (0.25)  (0.35)  (0.505)  =  0.50 


and  the  dimensional  2-percent  runup  is 

Ru2%  =  0.50  h  =  0.50  (20  ft)  =  10.0  ft 


The  present  Coastal  Engineering  Manual  method  given  by  Equation  4  estimates  runup  to  be 
Ru2%=  11.5ft. 


SUM  MARY :  This  CHETN  has  described  a  new  empirical  formula  for  estimating  the  vertical  runup 
distance  above  the  swl  that  will  be  exceeded  by  only  2  percent  of  the  irregular  wave  runups  on 
rough,  impermeable  slopes.  The  formula  is  based  on  the  same  hypothesis  used  in  an  earlier  analysis 
of  runup  on  smooth,  impermeable  slopes  (Hughes  2003b,  2004b),  i.e.,  the  weight  of  water  above  swl 
at  maximum  runup  is  proportional  to  the  maximum  depth-integrated  wave  momentum  flux  occurring 
in  a  wave  just  before  it  reaches  the  toe  of  the  impermeable  plane  slope. 

The  irregular  wave  runup  data  on  rough/riprap  slopes  of  Ahrens  and  Heimbaugh  (1988)  and  van  der 
Meer  and  Stam  (1992)  were  analyzed  in  terms  of  the  nondimensional  wave  momentum  flux 
parameter.  There  was  only  slight  difference  in  runup  between  waves  that  broke  on  the  slope 
(plunging  and  spilling  breakers)  and  nonbreaking  (surging  and  collapsing)  waves.  The  laboratory 
data  were  well  represented  by  applying  a  constant  reduction  factor  to  the  corresponding  equation  for 
runup  on  smooth,  impermeable  slopes.  The  new  formula  (Equation  13)  has  reasonable  predictive 
capability  for  irregular  wave  runup  on  slopes  in  the  range  1:2  to  1:4.  An  example  calculation 
illustrates  application  of  the  runup  equation. 

POINTS  OF  CONTACT :  This  CHETN  is  a  product  of  the  Coastal  Inlet  Structures  work  unit  of 
the  Coastal  Inlets  Research  Program  (CIRP)  being  conducted  at  the  U.S.  Army  Engineer  Research 
and  Development  Center,  Coastal  and  Hydraulics  Laboratory.  Questions  about  this  technical  note 
can  be  addressed  to  Dr.  Steven  A.  Hughes  (Voice:  601-634-2026,  Fax:  601-634-3433,  email: 
Steven.A.Hughes@erdc.usace.army.mil).  For  information  about  the  Coastal  Inlets  Research  Program 
(CIRP),  please  contact  the  CIRP  Technical  Leader,  Dr.  Nicholas  C.  Kraus  at 
Nicholas.C.Kraus@erdc.usace.army.mil.  Beneficial  reviews  were  provided  by  Dr.  Jeff  Melby, 
Coastal  and  Hydraulics  Laboratory;  and  Mr.  John  Ahrens,  retired  Coastal  and  Hydraulics  Laboratory 
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and  Sea  Grant.  Special  thanks  to  Mr.  John  Ahrens  and  Dr.  Jentsje  van  der  Meer,  INFRAM,  for 
providing  their  original  irregular  wave  runup  data. 
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